Semiclassical Collapse of a Sphere of Dust and Hawking Radiation 
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The independent dynamical variables of a collapsing homogeneous sphere of dust are canonically quantised 
and coupled equations describing matter (dust) and gravitation are obtained. The conditions for the validity of 
the adiabatic (Born-Oppenheimer) and semiclassical approximations are derived. On neglecting back-reaction 
effects, it is also shown that in the vicinity of the horizon and inside the dust the Wightman function for a 
conformal scalar field coupled to a monopole emitter is thermal at the characteristic Hawking temperature. 



1. The semiclassical collapse 

The total classical action for a self gravitating 
sphere of homogeneous dust immersed in the vac- 
uum (Schwarzschild) space-time is taken to be a 
generalization of the Oppenheimer-Snyder model 
PP (c = K Boltzmann = 1, K = 8ttG = l 2 /fi): 
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where the dust is described by the scalar field 
4> = 4>{t) (to allow for change in the number of 
particles in the second quantized development), 
fi = m/h is the inverse of its Compton wave 
length and K c — Kod v h e (ri), with 
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is the conformal factor of the interior (Robertson- 
Walker) metric with positive, negative and zero 
curvature respectively, 
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The constant Kg is related to the mass parameter 
M of the external Schwarzschild metric by the 
junction condition M — p 3 K at the surface of 
the sphere r = p K c g . 



Canonical quantisation leads to the Wheeler- 
DcWitt equation for the total wave function ^ = 
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One may now choose a suitable operator ordering 
in the gravitational kinetic term and, following 
a previously employed procedure Q] (analogous 
to the Born-Oppenheimer approximation used in 
molecular dynamics), express $ in the factorized 
form 

*(K,<i>) = K$(K)x(4>,K), (5) 
where %j) and \ satisfy the following equations || : 
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and (A) = J d<p\* (K> <P) A x(K, tj>) for any oper- 
ator A. 

1.1. The dust wave function 

If we consider the semiclassical (WKB) approx- 
imation to the wave function tj), one has 
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where S e ff is the effective action satisfying the 
Hamilton-Jacobi equation associated with the 
L.H.S. of Eq. In such a semiclassical limit 

4> will be peaked at the classical trajectory K c . 
One may then define a (conformal) time variable 
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Further if the R.H.S. of Eq. (@) is small one gets 
from Eq. (Q) the Schrodinger equation for a har- 
monic oscillator 

K c H M Xs = ih^ , (10) 

where we have scaled the dynamical phase, Xs = 
X exp {- £ /" < H M (rf) ) K c dr,'}. In the adia- 
batic approximation for i^T c the solutions are S: 

X, Y =7V JvJ ff JV ( 7 ^)e-^^/ 2 , (11) 

where iJ N is the Hermite polynomial of degree 
N, N N a normalization factor and 7 2 = /i K c 3 /h. 
The expectation value of the matter Hamiltonian 
on the solutions is given by 

K c {H M ) N = {N + l/2)h l iK c . (12) 

One can build a coherent superposition of the 
states in Eq. ( pi] ) corresponding to a minimum 
wave packet oscillating about tj> = with clas- 
sical frequency h/j,K c , which is constant in the 
adiabatic approximation, thus leading to a clas- 
sical trajectory for the dust in the classical limit 
(&-0)||. 

1.2. The gravitational wave function 

It is straightforward to derive an expression for 
$ which satisfies Eq. (0) in the semiclassical limit 
with vanishing R.H.S. (fluctuations). Suppose 
that at time rj for which d n K c = the gravita- 
tional wave function is a gaussian packet centred 
at the value K c (r) ) with width b. At a succeeding 
time r], tp will be given by 

$(K,n) = J dK' G e (K, tj; K', r) Q ) ip(K'; tj ) ,(13) 

where G e is the Green's function for the (inverted) 
harmonic oscillator for e = +1(— 1) or for a par- 
ticle moving in a linear potential for e = 0. One 
then obtains [01 
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We further note that the total gravitational wave 
function K ip(K) is for K = 0, thus giving no 
finite probability to fall exactly into the point like 
singularity K c = 0. 

If one now considers the limit fi — > followed 
by b — * (classical point-like limit) the classical 
trajectory is obtained: 
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where ( ^ | _ftT 2 | ip ) is the norm of the complete 
wave function Kip. With the above limits inter- 
changed, that is & — > with h finite, one would 
obtain 
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However, we expect that when b becomes smaller 
than the Planck length l v quantum gravitational 
effects (fluctuations) become significant. Hence 
it is more sensible to consider the limit for which 
one has b ~ £ p (corresponding to a minimum size 
wave packet of the order of the Planck length) 
and then consider £ p — > 0, which again leads to 
Eq. (U). 

From Eq. ( |l2] ) and taking into account the junc- 
tion condition at the surface po one may conclude 
that the total mass M G = M/k of the collaps- 
ing sphere is proportional to the number of dust 
quanta TV, 
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The above apparently differs from the result ob- 
tained previously, M G ~ m p ^/n 0, but we ob- 
serve that the N appearing in Eq. ( |l8| ) is an en- 
ergy quantum number of the dust, while the n 
in Bekenstein's formula can be obtained in the 
present context from quantising the gravitational 
Hamiltonian with purely classical matter source 
in the e = +1 case only Q. 

1.3. Consistency conditions 

Even if ip(K, if) is peaked on the classical value 
K r , there are fluctuations around it and one must 
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require Ehrenfest's theorem to hold in order to 
allow for the junction condition. In our case this 
means that one must have 
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For the wave function ip in Eq. (|lj) one finds 
A ~ 1/a 2 — ► in the classical point-like limit, 
Eq. (|16|). However if one considers the limiting 
procedure used in Eq. jl7| ) one has A ~ 0.18. 

It is particularly interesting to consider small 
black holes. Let us then take K c small at fixed 
h and b ~ t v , one then has A ~ £ p /K c , which 
implies that one must have r (ry) = p K c (rj) 3> 
£ p , or, in the limit for which r approaches r H = 
2 M from outside, 



M > £ p ^ M G > m p 
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This coincides with what one would expect from 
a naive quantum mechanical argument: if kTi/M 
is the Compton wave length of the black hole, it 
should be much less than its Schwarzschild radius 
2 M for the semiclassical approximation to work. 

The remaining consistency conditions are that 
the R.H.S.s of Eqs. (||) and (^) be much smaller 
than ( H M ) (for fluctuations in the number of 
dust particles to be negligible) and that in Eq. (|7]) 
the second term on the L.H.S. be larger than the 
third (for the adiabatic approximation to hold). 
For N ^> 1, the above three conditions are es- 
sentially the same and, on using Eq. (fsb , one 
obtains the following result M 
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which means that the Schwarzschild radius of the 
dust must be much larger than the Compton wave 
length of the dust particles. Since m <C m p for 
all known elementary particles, the latter expres- 
sion is a stronger condition than Eq. ( po|) so that, 
when it holds, one can safely match the interior 
with an external classical space-time. The fact 
that the adiabatic approximation is not valid if 
Eq. ( pl| ) is not satisfied is not surprising. Indeed 
in such a case fluctuations, corresponding to the 
creation of matter particles, are large Q (and the 
evaporation time small ||). 



2. Hawking radiation 

Let us consider the outer dust shell of the 
sphere situated at p Q . From the point of view 
of a static distant observer the position of the 
shell is not a classical variable but is a quantum 
observable determined by the gravitational wave 
function describing the semiclassical collapse ob- 
tained in the previous section, Eq. (pT[). 

We consider an isotropic massless scalar field 
ip = (p(p, rf) conformally coupled to gravity |t]] 
and to a static emitter. Its Lagrangian density 
will be given by 
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where R is the curvature scalar, Q{p, r) describes 
a particle (monopolc) emitter and the factor 5{p— 
p a ) \Ktp\ 2 /(ip \ K 2 | -0) forces the interaction to 
be localized on the outer shell. 

Let us assume the emitter has a discrete set 
of internal energy eigenstates described by | E ) 
where E is the energy observed by our distant 
observer. We suppose the emitter is initially in an 
excited state E and decays by emitting quanta 
of the scalar field (p to a state E — Kui. One 
may then estimate, using first order perturbation 
theory, the total probability amplitude P(u>,fj) 
for the emitter to decay in a finite conformal time. 
It will be given by 
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where Q w = \{E | Q(p o ,0) | E - hu)\ is the 
absolute value of the matrix element of the 
monopole between the states with energy Eq and 
E — Tluj, t ~ — 2 M \n[r}—T} H ) is the Schwarzschild 
time measured by the static distant observer ex- 
pressed in term of the conformal time, r(r} B ) = r H 
and fj is an upper cut off such that r\ H < fj < 
T] . Further during the interval (% , fj) the static 
monopole emitter near the horizon is immersed in 
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the gravitational wave packet associated with the 
last dust shell and D+ is the Wightman func- 
tion for the isotropic conformal scalar field in 
the Robertson- Walker metric ]tJ evaluated at the 
same spatial point p" = p' = p, but with the K c 
dependence "smeared" by the wave function Kip. 

One obtains a probability amplitude per unit 
(Schwarzschild) time equal to 0] 
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which is a Planck distribution with the usual 
Hawking temperature T H = 1/(3 = h/8irM and 
R(po) is a non zero constant given by R(r) H ) with 
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(see Ref. Q for explicit expressions). 

In the classical point-like limit, Eq. (fl6|), R —* 
1 and one recovers the usual field theory in the 
fixed Robertson- Walker-like background (H|7). 
However, in alternative limit, Eq. (|l7|), one finds 
that the emitter decouples from the conformal 
field, since R ~ a — ► 0, and there is no emission 
within any finite time T, 



P(w,T)-0. 



(26) 



One may speculate that this effect can be used 
to eliminate ultra-planckian effects. In fact, it 
is known that, if u> is the frequency of the emit- 
ted quanta as is measured by a distant observer, 
a fixed observer located near the point of emis- 
sion at r = r will measure instead a blue-shifted 
frequency uj* — (1 — 2 Af/r ) -1 / 2 to which gives 
uj* > m p /h — Ip^ 1 for r sufficiently close to 
r H . In order to create these modes, one must 
use an emitter localized in a region smaller than 
lu^ 1 ~ £ p . So one expects that b in the wave 
function Eq. ( [Til) should be less than l v for our 
collapsing shell to couple with conformal quanta 
of ultra-planckian energies and this would cor- 
respond to the limit in Eq. ( [l7|) which in turn 
implies Eq. (|26|). 

A further point worth noting is that R(oo) ac- 
tually depends on e through a (see Eq. (jig)). 
This of course implies that the probability am- 
plitude in Eq. (p3) depends on the internal dust 



sphere metric and one has different emission rates 
depending on the internal geometry. 

Let us end with a speculation: we consid- 
ered the matching condition in the classical limit; 
if one wished to consider quantum mechanical 
corrections one should replace it by the conti- 
nuity of the gravitational wave function inside 
the dust (Robertson- Walker) with that outside 
(Schwarzschild). This suggests the results in 
Eq. ( p4| ) should not change dramatically immedi- 
ately outside the dust, implying a form of quan- 
tum hair leading to information on the geometry 
in the black hole through the intensity of the ra- 
diation. 
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